Weak Mott insulators on the triangular lattice: possibility of a gapless nematic 
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We study the energetics of Gutzwiller projected BCS states of various symmetries for the tri- 
angular lattice antiferromagnet with a four particle ring exchange using variational Monte Carlo 
methods. In a range of parameters the energetically favored state is found to be a projected d x z_ y 2 
paired state which breaks lattice rotational symmetry. We show that the properties of this nematic 
or orientationally ordered paired spin liquid state as a function of temperature and pressure can 
account for many of the experiments on organic materials. 
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In the last few years the quasi-two dimensional organic 
salts k - (ET) 2 Cu 2 (CN) 3 and EtMe z Sb[Pd{dmit) 2 ] 2 
(abbreviated respectively as nCN and DM IT in the pa- 
per) have emerged as possible realizations of Mott in- 
sulators in the longsought 'quantum spin liquid' state 
0, 0, 0, 0, 0, 0, 0,11 • These layered materials are be- 
lieved to be well described by the single band Hubbard 
model on a nearly isotropic triangular lattice. At ambi- 
ent pressure they are Mott insulators which do not or- 
der magnetically down to temperatures ~ 30mA' (much 
lower than the exchange J ~ 250A inferred from high 
temperature susceptibility) 0, 0] ■ The low temperature 
phase is characterized by a linear T dependent heat ca- 
pacity, and a finite spin susceptibility just like in a metal 
(even though the material is insulating) [l], 0, 0] indicat- 
ing the presence of low lying spin excitations. There is 
a sharp crossover or possibly a phase transition at a low 
temperature « 5 A signalled by a peak in the heat capac- 
ity, and the onset of a drop in the susceptibility [l|, [j] . 
Further an external magnetic field induces inhomogenc- 
ity that is evidenced by a broadening of the NMR line d . 
Application of moderate pressure (« 0.5GPa) induces a 
transition to a superconductor (kCN) or metal (DM IT) 



At this point several questions arise: What is a good 
description of the spin liquid Mott state? What is 
the connection between the superconducting state and 
the underlying spin liquid state that becomes unstable 
upon applying pressure? What is the nature of the fi- 
nite temperature transitions/crossovers? We find, using 
a variational Monte Carlo analysis of the energetics of 
several possible wave functions for a spin Hamiltonian 
with Heisenberg and ring exchange interactions, that the 
nodal d-wave projected BCS state is the best candidate 
for the spin liquid. This state has gapless spin excita- 
tions and can naturally explain many of the experiments 
in kCN though a number of open questions remain. 



Our results are based on the model Hamiltonian d 

H = 2J 2 S r .Sr' + J4^( P 1234 + /l.c) (1) 

<rr'> □ 

Here S r are spin- 1/2 operators at the sites of a triangular 
lattice. The second term sums over all elementary par- 
allelograms, and P1234 performs a cyclic exchange of the 
four spins at the sites of the parallelogram. The multi- 
ple ring exchange is expected to be significant due to the 
proximity to the Mott transition in the organics. It is 
known that the 3 sublattice Neel order vanishes beyond 
a critical J4/J2 ~ 0.1 [l(| that can lead to novel spin 
liquid phases with no long range spin order. 

Previous studies on the above model, used a Gutzwiller 
rejected filled Fermi sea to interpret the experimcntsl 
The low energy theory of this state is described 
by a gapless Fermi surface of neutral spin-1/2 fermionic 
spinons coupled to a massless U(l) gauge field (also ob- 
tained [ll| within a Hubbard model description). The 
sharp crossover at T ~ 5K was associated with 'pairing' 
of spinons (l2l | . The condensation of the spinon pair field 
gaps out the t/(l) gauge field, and the resulting state is 
described as a Z 2 spin liquid. Ref. [13 suggested an exotic 
paired state that retains a finite gapless portion of the 
spinon Fermi surface. The possibility of a more conven- 
tional triplet paired Z 2 state induced by Kohn-Luttinger 
effects UJ of the spinons has also been pointed out. 

Here we study various paired spin liquid states for the 
J 2 — J4 model using variational Monte Carlo calculations. 
In terms of wavefunctions, paired states may be described 
by Gutzwiller projected BCS states. Two natural states 
(which retain the full symmetry of the triangular lat- 
tice) are projected singlet d x i_ y ? +id xy and nodal triplet 



2 wave states. Remarkably in a range of J4/J2 

state has bet- 
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we find that a projected singlet d x 
ter energy than either of these states. The d x 2_ y 2 state 
is a gapless Z 2 spin liquid state with nodal fermionic 
spinons and gapped Z 2 vortices (visons). In addition it 
spontaneously breaks the discrete rotational symmetry of 
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the triangular lattice but preserves lattice translational 
symmetry. Thus it is a gapless Zi spin liquid coexisting 
with a 'nematic' or orientational order parameter. The 
pairing structure of the spinons determines the pairing 
structure of the superconductor that forms under pres- 
sure. Thus we propose a nodal ci-wave state for the pres- 
sure induced superconductor as well. Due to the discrete 
broken rotational symmetry both the insulator and the 
superconductor will have non-trivial finite temperature 
phase transitions in an ideal sample. We describe these 
and comment on their implications for the experiments. 

Various variational states may be constructed by start- 
ing with a system of spin-1/2 fcrmionic spinons f ra hop- 
ping on a finite triangular lattice of size LI x L2 at half- 
filling with a "mean field" Hamiltonian: 



Hmf = 



E 



-trr-fUr'c + (Kt> f^ft- 1 + h.C.)] (2) 



The variational spin wave-function |^)„ ar = Pg\^)mf 
where the Gutzwiller projector Pq = FT. (1 — rijtnji) 
ensures exactly one spinon per site. Unknown pa- 
rameters in l^uar are fixed by minimizing the energy 

E V ar = (*«or|-H"|*)'Uor/(*«ar|*)Kor (with H given by 

Eq. [1]) with only nearest neighbor t rr > = t. The 
simplest \*&)mf corresponds to A rr > = i.e. a filled 
Fermi sea. The corresponding \^) va r = \^)pfl = 
Pa {jXka /| ) |0) with no variational parameters. The 

prime on the product implies restriction to k such that 
the singlc-spinon level eg < Ef, the Fermi energy. More 
complex variational wave-functions are obtained with 
different patterns of non-zero A rr i which correspond 
to various projected BCS wavefunctions \^}pbcs = 

N/2 

|0). Here <fe = 



P G \BCS) = Pa^hflf^ 
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with £g = eg — /i. Further we 

write Ag, the Fourier transform of A rr /, as Ag = Ao-F(fc) 
where the form of F(k) is fully determined from a partic- 
ular pattern of A rr i (or equivalently, a particular Cooper 
pair channel). The two variational parameters: gap pa- 
rameter Ao and the 'chemical potential' /x, are both de- 
termined by minimizing the energy. 

The properties of the three superconducting gap func- 
tions are as follows: The d x 2_ y 2 + id xy state is invariant 
under spin rotation, lattice rotation and translation sym- 
metries, but breaks both time-reversal and parity. Af- 
ter projection it corresponds to the 'chiral spin-liquid' 
state [14J. The d x n_ y n state is a spin-singlet with a 
cos(29) angular dependence (8 is the angle subtended 
by a bond). It breaks lattice rotational symmetry while 
preserving translations and time-reversal. Finally, the 
triplet / x 3_ 3l!/ 2-wave state has orbital part varying as 
cos(39) while in spin-space it has zero projection along 
a quantization axis. This breaks spin rotation but pre- 
serves all the lattice symmetries and time reversal. Both 
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FIG. 1: (a) Difference between Epbcs and Epfl for various 
paired states in the units of J2. The calculations are done on 
a lOx 11 lattice with anti-periodic boundary conditions using 
the standard Metropolis Monte Carlo [3 with 10 5 sweeps, 
(b) Gap parameter Ao for various paired states in units of t. 



and f x 



-3xy 



2 possess nodes along the Fermi sur- 



face in fc-space while d x 2_ y 2 + id xy is fully gapped. 

Figure QTa) shows the the difference 5E = Epbcs — 
Epfl for the three paired states, namely projected 
d x 2 —y 2 ~t~ id x y-} fx 3 -3xy 2 an d d x 2_ y 2. Clearly for J4/J2 ^ 
0.25 the projected Fermi liquid is the best variational 
state. Interestingly, for a wide range of parameters 
0.10 < J4/J2 ^ 0.23 the projected d x 2_ y 2 wins over all 
projected Fermi liquid as well other two paired states. In 
the regime 0.23 < J4/J2 ^ 0.25, the error bars preclude 
any conclusion. Figure QIb) shows the optimal value of 
the gap parameter A for these three states. Consis- 
tent with the results for optimal energy, Ao ~ for 
J4/J2 > 0.25 while for 0.10 < J4/J2 < 0.23 the state 
d x 2_ y 2 has a non-zero and largest value of Ao among all 
paired states. For smaller values of Ja/ 'hi^, 0.10) it is 
expected that the spin-rotation symmetry breaking spiral 
state would be the ground state of H [lfj. In addition, 
we also studied a projected Fermi liquid with staggered 
flux $ through alternate triangular plaqucttes. We found 
that the energy has minima at $ = 0, ir and that the tt 
flux state is always higher in energy than the zero flux 
state \^)pfl for all values of J4/J2. 

Heuristically, large values of J4/J2 favors delocaliza- 
tion of electrons. Thus it is not surprising that PFL 
is the ground state for large J±/J2- Since J2 > 0, the 
triplet paired /-wave state is expected to be unfavor- 
able, consistent with our results. Further, the electrons 
are more delocalized in the uoAdl-d x 2 _ y 2 compared to 
d x 2_ y 2 + id xy since the latter is fully gapped. Thus for 
values of J4/J2 not so small as to induce spiral order for 
spins but small enough that PFL is destabilized, our re- 
sult that a projected nodal paired state is favored seems 
reasonable. Our results connect well with earlier varia- 
tional Monte Carlo 2, 16| and other numerical studies 



1 71 ] of superconducting states in anisotropic triangular 
lattice Hubbard models which also found good evidence 
for a nodal d-wave state. 
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FIG. 2: Schematic sketch of the expected specific heat curve. 
The left panel shows the "mean field" behavior. The full line 
is in the clean limit, and the dashed line the modification due 
to impurity scattering. The right panel shows the expected 
behavior when gauge fluctuations beyond mean field are in- 
cluded. 



We now describe various properties of the state de- 
scribed by the nodal d-wave paired wavefunction. A 
mean field Hamiltonian that describes the excitations of 
this state is simply the Hmf of Eqn. [2 Fluctuations 
about the mean field state are described by coupling the 
spinons to a Z 2 gauge field. This state is thus an exam- 
ple of a Zi spin liquid. The excitations of the Z 2 gauge 
field are Z 2 flux configurations (known as visons) which 
are gapped in this spin liquid phase. The low energy 
physics is then correctly described by the BCS Hamil- 
tonian Hmf- Many properties of the nodal d-wave spin 
liquid at low temperature are thus similar to the familiar 
spin physics of a nodal d-wave superconductor. We now 
describe some of these in relation to the experiments. 

Specific heat and spin susceptibility: In the absence of 
impurities, the density of states for a nodal supcrfluid 
vanishes linearly with energy and consequently the spe- 
cific heat C = aT 2 where the coefficient a is, in principle, 
determined by the velocities that characterize the nodal 
dispersion of the spinons. Impurity scattering generates 
a non-zero density of states leading to a specific heat 
C ~ jT, and a constant spin susceptibility Xo as T — > 0. 
Further the low-T Wilson ratio xT jC is constant of order 
one. All of these are in agreement with the experiments 
on the organic spin liquid materials. 

The impurity scattering rate can be roughly estimated 
by equating the entropies of the paired nodal spin liq- 
uid and 'normal' states at T* , where T* is a mean field 
or crossover scale below which the pairing sets in (Fig. 
[5]) . Above T* spinons may be described as having a gap- 
less Fermi surface with a specific heat C = jspinonT. 
Here ^/spinon = (jf /3)fe^^spmon(-^/) and Ti S pi non (^Ef^ = 
0.28 /t S pi non is the spinon density of states at the Fermi 
energy. Equating the entropies of the paired and nor- 
mal states at T*, we estimate a = 2 r y sp i non /T* . Impu- 
rities will cut-off the T 2 specific heat of the nodal spin 
liquid at a scale T and lead to a low-T gamma coeffi- 
cient given by 7 — aT/k B = 2^ spinon T / (k B T*). Now 
tspinon ^5 2J 2 ~ 250 K and from the measured low-T 



specific heat with 7 « 15m JK ~ 2 moZ _1 , we estimate an 
impurity scattering rate Y < 0.25k B T* w 1.5K. While T 
is reasonably small compared to superconducting gap it 
is appreciable enough to generate a constant density of 
states at low energy and lead to an apparent Fermi-liquid 
like behavior in the specific heat at the lowest tempera- 
tures accessible, consistent with the experiments 0] ■ The 
experiments also apparently show that the low tempera- 
ture linear specific heat is insensitive to magnetic fields 
upto about 8 T. This poses a difficulty for the present 
theory as the Zeeman coupling to the field is expected 
to increase the low energy density of states. Indeed the 
prior proposal of an 'Amperean' paired state [l2T | was 
partly motivated by the insensitivity of the specific heat 
to a magnetic field. However as discussed below the Am- 
perean pairing has some difficulty with describing the 
superconducting state that develops under pressure. 

Thermal conductivity: Nodal spinons (as also dirty d- 
wave superconductors) lead to a finite 'universal' metallic 
thermal conductivity k ~ T as T — > 0. In practice how- 
ever observation of this effect requires low temperatures 
to eliminate the phonon contribution. In thermal trans- 
port measurements on kCN, a plot of k/T as a function 
of T 2 of data above 1 K indeed extrapolates to a constant 
in the zero temperature limit. However, data for T < 0.5 
K rapidly extrapolates to zero and has been interpreted 
as evidence for a gap [l9| . We do not have an explanation 
of this phenomenon. 

Field induced inhomogeneity: At ambient pressure 
NMR studies of kCN show the development of a 
magnetic-field induced inhomogeneity Within our 
theory this may be rationalized as follows. Due to the 
effect proposed by Motrunich [2(j , the external magnetic 
field induces an internal magnetic field for the spinons 
which can lead to vortices (visons) of the spinon pair con- 
densate. The resulting 'mixed state' is inhomogeneous 
that occurs below the pairing scale and increases in pro- 
portion to the field. 

T = phase diagram under pressure and superconductivity: 
In general pressure increases the ratio t/U which implies 
an increase of J4/ J 2 leading to suppression of the pair 
amplitude. Thus increasing pressure suppresses the 
pairing transition. Increasing pressure also leads to an 
insulator-to-metal transition. Clearly two situations 
are possible depending on whether the pair order is 
killed before or after this metal-insulator transition. 
In the latter case superconductivity will be obtained 
in the metal close to the Mott phase boundary. We 
propose that this is realized in kCN. On the other hand 
superconductivity has not been found in DMIT under 
pressure. We suggest that in this material the pair 
order is killed under pressure before the metal-insulator 
transition. An interesting experimental test of this 
suggestion is to study the Mott insulating phase of 
DMIT at pressures just below the metal-insulator phase 
boundary. Here the spinon Fermi surface state, with its 
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characteristic signatures such as, for instance, the T 2 / 3 
heat capacity (produced by gauge fluctuations) will then 
survive to low-T without any pairing transition. 

If the pairing extends into the metallic phase the su- 
perconductor that results will also have d x 2_ y 2 symmetry, 
and will (for an ideal isotropic triangular lattice) break 
lattice rotational symmetry. The spinons of the insulator 
now become the nodal Bogoliubov quasiparticles of this 
d-wave superconductor. Thus the low temperature spe- 
cific heat and spin susceptibilities of the superconductor 
will behave similarly to that of the spin liquid insula- 
tor. Further the NMR relaxation rate 1/T X T ~ T 2 for 
T > r (the impurity scattering rate) and will saturate 
to a constant at the lowest temperatures. The former is 
in agreement with existing data on kCN for T close to 
T c [5j . Such a relaxation rate is not expected within the 
alternate Amperean paired state making it difficult 
to connect the pairing transition in the spin liquid with 
that in the metal. The NMR data [H also shows that 
the Knight shift is only weakly suppressed on entering 
the superconducting state. However this may be due to 
complications associated with sample heating 21 1 . 
Finite-T phase diagram: The d x i_ y i paired spin liquid 
breaks a discrete orientational symmetry. As the temper- 
ature is increased, it undergoes a phase transition in the 
three-state Potts universality class at a certain temper- 
ature. The nematic order also leads to a richer finite-T 
phase diagram in the superconductor. A straightforward 
analysis shows that at any non-zero T, the superconduct- 
ing order becomes power law while the discrete broken 
orientational symmetry is preserved up to a certain tem- 
perature, say, T\ after which the system enters a 'float- 
ing' phase where both supcrfluid and orientational orders 
have power law correlations. Finally at an even higher 
temperature, say, T2 there is a further phase transition to 
the normal state with exponentially falling correlations 
for both orders. Both the phase transitions at T\ and T2 
are of the Kosterlitz-Thouless type with very weak signa- 
tures in the specific heat. However in both the insulator 
and the superconductor for the initial pairing transtion, 
by the usual Ginzburg criterion, the fluctuation regime 
will be rather small. In practice, there is a small lattice 
anisotropy in both kCN and DM IT which will pin the 
nematic order parameter and smoothen out any sharp fi- 
nite temperature nematic transition. Further even in an 
isotropic material weak disorder acts as a 'random field' 
on the nematic order parameter and will kill the finite 
temperature nematic phase transition. Nevertheless if 
the disorder and the lattice anisotropy are weak, a sharp 
crossover behavior associated with the paired nematic or- 
der might be expected in the insulator. Such a crossover 
is visible in the existing experiments in a variety of prop- 
erties. 

Thus our proposal of a nodal d-wave Z2 spin liquid 
state that develops out of the spinon Fermi surface spin 
liquid can account for a number of measured properties 



of the organic materials. The most serious difficulty at 
present is the field independence of the specific heat. Per- 
haps the impurity scattering rate T can be bigger than 
our rough estimate without any major effect on the pair- 
ing transition temperature. 

Are there any 'smoking gun' tests of our proposal for 
future experiments? The essence of our proposed gap- 
less Z2 spin liquid state can be probed in experiments 
on kCN through the flux trapping effect described in 
Ref. . We note that such experiments have been per- 
formed on the cuprate materials [23| , and may be feasible 
in the organics as well. A Josephson tunneling exper- 
iment with a spin liquid insulating barrier between two 
superconductors has also been proposed as a probe of the 
gapped charge-e spin-O charge carriers expected in this 
spin liquid state 24 1. 

To summarize we studied the energetics of various 
Gutzwiller projected BCS states for the triangular lat- 
tice antiferromagnet with a four particle ring exchange. 
In a range of parameters the best state is a projected 
d x 2_ y 2 paired state which breaks lattice rotational sym- 
metry We described many properties of this state in 
light of experiments on the candidate organic spin liquid 
materials. Future experiments will hopefully shed light 
on whether such a paired nematic spin liquid really exists 
in these materials. 
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